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Abstract 

We present a new closed form for the interpolating polynomial of the 
general univariate Hermite interpolation that requires only calculation of 
polynomial derivatives, instead of derivatives of rational functions. This 
result is used to obtain a new simultaneous polynomial division by a com- 
mon divisor over a perfect field. The above findings are utilized to obtain 
a closed formula for the semi-simple part of the Jordan decomposition 
of a matrix. Finally, a number of new identities involving polynomial 
derivatives are obtained, based on the proposed simultaneous polynomial 
division. The proposed explicit formula for the semi-simple part has been 
implemented using the Matlab programming environment. 

1 Introduction 

The Hermite interpolation of total degree is described in the following Theorem 

Theorem 1.1. Given n distinct elements Aq, Ai, . . . , A„_i in a perfect field k, 
positive integers rui, i — 0, . . . ,n — 1, and aij G k for 0<i<7i — 1, 0<j< 
mj — 1, then there exists one and only one polynomial r G k[a;] of degree less 
than X]r=o '^i> such that 

r(j)(Ai) = fly, < j < m, - 1, 0<i<n-l. (1.1) 
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This polynomial r is explicitly given by, 



n— 1 rrii — l m^— j — 1 



Ux-X,) 

ji k\ [ n{x) 



r{x) 



i=0 j=0 fc=0 



n{x) 



{x - XiY"--!-'"' 




n{x) - n(x- AO 



and r('^'(a) is i/ie k—th derivative of r at a. 



The applications of hermite interpolation to numerical analysis is well known. 
In this paper we show the use of a computationally more efficient formula for 
the Hermite interpolation to a number of algebraic applications. 

A number of forms of the interpolating polynomial r{x) have been reported 
in the literature, however they all require calculation of derivatives of rational 
polynomial functions (e.g. [T]), or recursive calculation of the coefficients aij of 
Theorem 11.11 (e.g. [S]). We firstly propose a new closed form of the interpo- 
lating polynomial r of the univariate Hermite interpolation that requires only 
calculation of polynomial derivatives instead of derivatives of rational functions, 
as required in Theorem 11.11 thus achieving a substantial computational accel- 
eration compared to the Theorem ll.il without imposing any restrictions in the 
generality of the interpolation. We then exploit this result to show a Theorem 
for the calculation of the remainders of the simultaneous division of a number 
of polynomials by a separable polynomial divisor, as derivatives of a unique 
polynomial. Namely, let us note that the remainder r € lk[a;] of the Euclidean 
division of any polynomial P G k[x] of degree n by a separable polynomial 
Q S k[a;] of degree m, where n > m, can be calculated in closed form using the 
Langrange interpolation formula as following: 



where Ai , . . . , Am are the roots of Q in the algebraic closure Ik of k. In this 
work we extent this simplified idea of polynomial remainder calculation by Lan- 
grange interpolation, to achieve simultaneous polynomial division by a common 
separable divisor, using the closed form of the interpolating polynomial of the 
univariate Hermite interpolation r that we propose. 

Furthermore, the above results will be used to obtain a new closed formula 
for the semi-simple part of the Jordan decomposition of an algebraic element 
in an arbitrary algebra. If k is a perfect field and A an algebra over k with unit 
1, then the well-known result on the Jordan decomposition, (e.g. [2l|4l[7]) is 
presented in the following theorem. 
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Theorem 1.2. For each algebraic ^ G A there exist unique Sa,Na G k[A] 
such that A = Sa + Na, Sa is semi-simple and Na is nilpotent. 

It should be noted that in the case of k being algebraicahy closed and A is 
a subalgebra of the matrix-algebra k"^", then Sa is semi-simple if and only if 
iS*^ is diagonalizable. 

Theorem 1 1.2 1 has various interesting applications, such as efficient computa- 
tion of high powers of A, therefore, efficient algorithms for the calculation of 
Sa in terms of A are of great importance. A proof of the existence of Sa that 
is presented in the book of Hoffman and Kunze [1], is based on the Newton's 
method and yields direct methods for computations. An algorithm, which is 
essentially based on these ideas, is given by Levelt [B]. The algorithm of Bour- 
goyne and Cushman [3] is faster, because higher derivatives are used. In [9j D. 
Schmidt has used the Newton's method to construct the semi-simple part of 
the Jordan decomposition of an algebraic element in an arbitrary algebra, show- 
ing quadratic convergence of the algorithm. Another approach uses the partial 
fractions decomposition of the reciprocal of the minimal polynomial [3 110( lllj . 
These works however, do not provide a closed form formula for Sa as a polyno- 
mial of A. An explicit construction of the spectral decomposition of a matrix 
using Hermite interpolation appears in |11| . However, this approach requires 
the use of Taylor coefficients of the reciprocal of the matrix minimal polynomial. 

In this work, we exploit our proposed simultaneous polynomial division by a 
common separable divisor, to obtain a new closed formula for the semi-simple 
part Sa of the Jordan decomposition of an algebraic element A in an arbitrary 
algebra. The proposed closed formula requires only evaluation of the derivatives 
of the basic Hermite-like interpolation polynomials that are associated by the 
eigenvalues of A, up to the maximum algebraic multiplicity of the roots of 
the minimal polynomial of A, as well as matrix multiplication operations. The 
derived closed formula for the semi-simple part of the Jordan decomposition has 
been implemented in the Matlab programming environment. The sourse code 
is provided and results are shown from it's executions on examples of matrices 
taken from the literature. 

Furthermore, a number of new identities involving polynomial derivatives are 
also shown, based on the the simultaneous polynomial division by a common 
separable divisor. 

The rest of the paper is organized as follows. In section 2 we present a 
new closed form for Hermite interpolation. In section 3 we show a Theorem 
for the calculation of the remainders of the simultaneous division of a number 
of polynomials by a separable polynomial divisor, as derivatives of a unique 
polynomial. In section 4 we present an explicit Formula for Sa, the relative 
source code in MATLAB environment, as well as some numerical examples from 
the literature. Finally, in section 5 the main result of section 3 is generalized 
and applied to the derivation of some interesting polynomial formulas. 
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2 A new closed form for Hermite interpolation 

Let Ao, Ai, . . . , A„_i be distinct elements in a perfect field k and mo, mi, . . . , m„_i 
be positive integers. Let us denote by Lk the polynomials given by 

Tl — 1 / -i \ . 

Furthernore, we denote byA^, 0<fc<n— 1 the m^ x m^ lower triangular 
matrices [kj] G k"**'^"**', <i,j < ruk — 1, given by 



(;)(Lft)(^-j)(Aft) if 0<j<i<mfe-l 

if < i < j < mfe - 1 



where (ift)*-* ■'\^k) is the derivative of order [i — j) of the polynomial Lk{x) at 
A/j. Thus Afe has the following representation 

(o)ifc(Afe) ... 

Q(Lfe)«(A,) (l)i.(Afe) ... 



(-'■-i)(L,)(--i)(Afe) {^\-'){Lk)^^>'-^\Xk) . . . {ZT-^Lk{\k) 

(2.2) 

For our purpose the following technical Lemmas are required: 

Lemma 2.1. The matrices A/;, < A; < n — 1 are invertible with A^^^ = 
Silo" i^^k ~ ^kY, where is the ruk x m/j unit matrix. 

Proof. Clearly, for < fc < n — 1 holds Lk{Xk) = 1- Therefore, all matrices A/j, 
< fc < n — 1 are invertible and unit diagonal lower triangular. Therefore, it 
follows (Imk - Afe)"'= = and consequently A^ I]™o~^(^"fc ~ ^^Y = Im^- □ 

Using the Leibnitz's rule for derivatives, wc easily get the following Lemma: 

Lemma 2.2. For <i,s < ruk — 1 and < j,t < n — 1, the following holds: 



( _ \ \s I if t 1^ j 

r Lt{x)Y^ { If t = j andKs . (2.3) 

0{LjY'-'\Xj) if t = jands<i 



Our proposed form of Hermite interpolation can now be presented in the 
following Theorem. 

Theorem 2.3. Given n distinct elements Ao, Ai, . . . , A„_i in a perfect field k, 
positive integers m^, i = 0, . . . , n — 1 , and atj G k for < j < n — 1, < i < 
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rrij — 1, then there exists one and only one polynomial r £ k[a;] of degree less 
than X]"=ro^ such that 

r(')(Aj) = fly , < j < n - 1, < ? < - 1. 

This polynomial r is explicitly given by, 



(2.4) 



ji-i 



(2.5) 



j=0 fc=0 



where the matrices Xj, Aj are given by 



(mj-1)! ^jy^) 



and 



Aj = [ aoj Oi 



a. 



rrij-lj 



Proof. It can be observed that p.Sp can be equivalently written in the fohowing 
form; 

where 









Clj 


















_ Crrij-lj _ 




^nij — Ij 



,0 < J < n - 1. 



(2.6) 



Now, by calculating the derivative of order i of the polynomial r at Xj and using 



(|2.3I) in Lemma [2^ we obtain 

rW(A,) = ^c,4:)(L,)(-'=)(A,), 



fc=0 



(2.7) 



for allO<j<n — l,0<i<m— 1. Taking into consideration the definition 
of Aj in p.2p , the system of equations p.7p can be rewritten in the following 
matrix form 



r(A,) 










Clj 




-A, 




r(™j-i)(Aj) 







< j < n - 1 



(2.8) 
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By substitutting ((^ in ((^ . we get 









^'(A,) 










(^rrij — Ij 



Hence, we derive that r(')(Aj) = a^^ , < j < n — 1, < i < — 1. The second 
equahty of (|2.5p is obtained by Lemma [^TTl 

Moreover, it can be easily confirmed that any polynomial (x — XjY Lj(x), 
0<j<n — f,0<i< rrij — 1 has degree less than '^"Zq rrii, and since r is a 
k- linear combination of these polynomials, we conclude that the degree of r is 
less than J2^=o ■ '-' 



3 Simultaneous division of polynomials by sep- 
arable polynomial 

At this point we are ready to present the following generalization of Euclidean 
polynomial division by a separable divisor, based on Theorem 2.3. 

Theorem 3.1. Let g G k[a;] be a separable polynomial and Ao,...,A„_i be 
the roots of g in the algebraic closure k o/ k . Then for any polynomials 
/o, /i, . . . , fm-i € k[a::], there exists unique r G k[a;] of degree less than mn 
and unique polynomials qo, qi, g,„_i G k[a:] such that 



f^ 



' gQi 5 < i < m — 1 . 



This result is optimal, in the sense that if m > I and g is inseparable, then this 
result is not true. The polynomial r is given by 



where 



(m-l)! 



-Lj{x) 



/o(A,) /i(A,) . . . /™-i(A,) 
and Lj,Aj are respectively as in i2.1\) . 112. 2\) by mp — mi = ... = m„_i = m. 

0, . . . , A„_i) as vector space over 

(3.1) 



Proof. Let k be the dimension of the field 
k, that is 

fc= [k(Ao,...,A„_i),k] 
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then there exist Ti,...,Tfc_i m k(Ao, . . . , A„_i) such that {1,ti, Tk-i} is a 
basis of k(Ao, . . . , A„_i) as a vector space over k. 

Now, using Theorem [231 by mo = mi = ... = mn-i m, we have that there 
is unique polynomial r S k(Ao, . . . , A„_i)[x] given by (|2.5p for mo = mi = ... = 
m„_i = m having degree less than mn such that 

ir)'-'H^s) = M^j) (3.2) 

for all < i < 771 — 1 , < j < n — 1. Therefore there exist polynomials 
Qi G k(Ao, . . . , A„_i)[a:], < i < m — 1 such that 

^'^=f.+9%- (3.3) 

Moreover, by (|3.ip we have that the dimension of k( Aq , . . . , A„_ 1 ) [x] as free mod- 
ule over k[a;] is k and {1, ti, Tk-i} is a basis of k(Ao, . . . , A„_i)[a:;] over k[a;]. 
Therefore the polynomials r, qi G k(Ao, . . . , A„_i)[a;] can be uniquely written in 
the following form: 

n-l 

r = r + '^Tsrs, (3.4) 

and 

ft = ft + XI "^sfti, < i < m - 1, (3.5) 

s=l 

where r, r-^., qi, qsi £ k[a::] with degr, degr^ < mn. 

Setting (|3.4p and p.5p in p.3p and using that {1,ti, t^^i} is hnearly 
independent over k[a::], we get 

r(^'=,/.+5ft (3.6) 

for all < « < m — 1. 

Now from (|3.6p we clearly get that r satisfies the identities (|3.2p and since 
the polynomial of degree less than mn satisfying the identities p.2p is unique, 
we conclude that r = r. 

Finally, suppose that Theorem l3.1l is true for one polynomial g € k[a;] having 
a root A S k of multiplicity > 1. Then there exist a polynomial go in k[a;] such 
that 

g{x) = (x - Xfgo{x). 
Applying Theorem 13.11 for g and /o = /i = 1 we obtain 

r{x)^l + {x-\fgo{x)qo{x), (3.7) 

r^^\x)^l + {x-\fgo{x)q^{x) (3.8) 

for some r, qp, 91 G k[a;]. 

Differentiating p.7p . and setting x = A in the resulting identity, as well as, 
in p.Sp we respectively get 

r(i)(A) = and r^^\X) = 1. 

That is not true. □ 
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Corrolary 3.2. Let g be as in Theorem \3.1\ Then for any f € k[x], c E k— {0}, 
m G N with m > 2 there exists unique r G k[a;] of degree less than mn such that 
r ~ f mod g and r' = crg^^^^ . 

Proof. According to Theorem 13 . 1 1 there exists unique r G k[a;] of degree less than 
mn, such that 

f^r,cf = r' mod.g, ...,c™"V = r(™-i) mod 5. (3.9) 

Now muhiplying the i— th relation of (j3.9|) by c and afterwards substracting 
from the result the i + 1— th relation we get 

(r' — cr)^'-* = mod g, 

for all < i < m — 2 and since g is separable we direct get the conclusion. □ 

4 An explicit formula for Sa 



We will use Theorem 13.11 to give a closed formula for the semi-simple part of 
the Jordan decomposition of an algebraic number of an algebra A over a perfect 
field k. 

Lemma 4.1. Let g G k[a;] be a separable polynomial and Aq, . . . , A„_i are the 
roots of g in k. Let f G k[a;, y] be a polynomial of two variables and m be a 
positive integer. Then there exists unique r G k[a;] of degree less than mn such 
that 

f{x,y) = r{x + y) mod/, 

where I C k[a:, is the ideal generated from the polynomials g{x), y'". Further 
the polynomial r is given by 

n-l 

r(a;) = ^X,ATi^^. ek[x], (4.1) 
i=o 

where the matrices Xj, Aj are given by 



L,{x) ^L,ix) . . . ^Y^^iQ' L,{x) 

d 



= [ /(A„0) ^/(A„0) . . . ^/(A„0) 

and Lj,Aj are as in \2.1]) . 112. 2\) by mo = mi = ... = m„_i = m. 
Proof. The polynomial / can be rewritten in the following form 

il dy 

for some h G k[a:,y]. 



fix, y)^Yl ^)y' + y)' (4-2) 
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Furthermore, according to Theorem 13.11 there exists unique polynomial r G 
k[x] of degree less than mn, given by (I4.1[) . such that: 



^-T/(a;, 0) = r*^*^ (x) mod 5(2;), < i < m - 1, 



Combining (li?^ with we get 

m — 1 



mod /. 



(4.3) 



(4.4) 



Furthermore, by using the Taylor formula we have: 



E 



y' = r(x + y) mod 



(4.5) 



Substituting ()4.5p in (|4.4p . completes the required proof. 



□ 



Let p S k[a;] be the minimal polynomial of an algebraic element A of an 
algebra A over k with unit 1. Let Ai, « = 0, 1, . . . , n — 1 be the distinct roots 



of p in k, and let fc;, i = 1,2, 



1 be their respective multiplicities. We 



denote p{x) := Y[i=o ~ ^i) ^^"^ m(jj) := maxjfco, . . . , fcn-i}. 

Theorem 4.2. The semi-simple part Sa of the Jordan decomposition 0/ A is 
given by Sa = f{A), where r[x) is the polynomial 



(4.6) 



wher 



A. 



Lj{x) 
[ A, 



1! 



and Lj , Aj are as in h2.1\] 



2.2]) by niQ = mi 



(m-l)! 



-Lj{x) 



nin-i = m{p). 



Proof. Since Sa is the semi-simple part of A and Na is the nilpotent part of 
A, the minimal polynomials of Sa and Na are respectively p{x) and a:™'^-*. So 
we have p{Sa) = and N^'^^^ — 0. Now if we apply Lemma [4.11 by choosing 
/(x, y) = X, g{x) = p{x) and m = rn{p), and taking into account that /(x, 0) = 
X, and ^/(x, 0) = for 1 < i < m — 1 we have that for the polynomial r 



given by (|4.6p holds: 

X — r{x + y) mod / , (4.7) 

where / is the ideal generated from the polynomials p{x) and y"^^P\ Finally 

setting x = Sa and y = iV^ in we get Sa = r-(S'^ + iV^) = ?'(A). □ 
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4.1 An implementation of the explicit formula for Sa in 
the MATLAB programming environment and some arith- 
metic examples 

In this subsection we present a computational implementation in the MATLAB 
programming environment for the calculation of the semi-simple part Sa of the 
Jordan decomposition of A. The MATLAB source code of function semisimple 
is provided, with appropriate exlanatory comments and pointers to the relevant 
Equations. Instead of using the interpolating polynomial given in Theorem 14.21 
we implemented the equivalent formulation for the interpolating polynomial of 
Eq. dlH), by setting 

[coj,cij, . . . ,Cm(p)_ij]^ = Aj^[Xj,0, . . ., 0]"^. 

The function's input arguments are the matrix A whose semisimple part is 
required, as well as the number of decimal digits that will be taken into con- 
sideration for calculating the multiplicity of the eigenvalues of A. Integer m{p) 
in Theorem 14.21 which is equal to the maximum multiplicity of the roots of 
the minimal polynomial of A, is calculated as the maximum number of non 
zero elements of the rows and columns of the Jordan decomposition of A. The 
main function (semisimple) uses two more functions that are also provided, 
one for the calculation of the Lagrange polynomial coefficients and one for the 
calculation of the coefficients n*'' power of a given polynomial. 

function R=semisimple(A,tol) 

% by A. Kechriniotis, K. Delibasis, 2009 

% A: matrix whose semisimple part is required 

°/o tol: number of digits to use for finding multiple eigenvalues of A 

p=poly(A) ; 
a=roots (p) ; 

7o determine multiplicity of roots according to tolerance tol 

tol=10~tol; 

rl=round(a*tol) /tol ; 

[r_uniq, i_uniq] =unique(rl) ; 

for i=l : length(r_uniq) 

multi_uniq(i)=sum(rl==r_uniq(i) ) ; 

end 

a=r_uniq; 
n=length(a) ; 

°/o determine m as the max root multiplicity of the minimal polynomial 
J=jordan(A) ; 

m=max(max(sum( J~=0) ) , max(sum(J~=0,2))) ; 
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y, create a table with n rows to hold Langrange polynomial coefficients 
as in (2.1) 
for k=l:n 

myLangr_coef s(k, : )=myLangr (a,k,n) ; 

end 

"/o Calculate cnm array in (2.6) 
for jO=l:n 

generate Lamda mxm matrix, as in (2.2) 
for i=l:m 

for j=l:m 
if j>i 

Lamda(i, j)=0; 

else 

ll=myLangr_coef s(jO, : ) ; 

ll=polypower (11 ,m) ; 
f l=nchoosek(i-l , j-1) ; 

if i>l calc (i-1) order derivative of 11 

lk=ll; 

for k=l:i-j 

lk=polyder(lk) ; 

end 

Lamda(i, j)=f l*polyval(lk,a(jO)) ; 
else no derivative of 11 

Lamda(i, j)=f l*polyval(ll,a(jO)) ; 

end 

end 

end 

end 

I=diag(ones (size (A, 1) , 1) ) ; 
coll=[a(jO) ,zeros(l,m-l)] ' ; 
c(: , jO) = (Lainda"-l)*coll; 
coll ; 

end 

R=zeros (size (A) ) ; 
for i=l:m 

for j=l:n 

Ln=myLangr_coef s (j , : ) ; 

LnM=polypower (Ln,m) ; 

R=R+c(i, j)*polyvalm(LnM, A)* ((A-a(j)*I) ~ (i-1)) /factorial (i-1) ; 

end 

end 
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function Ln=myLangr (a,n,N) 
7o a: is the vector of roots 
7. N: length of a 

7„ n: integer between 1 and N inclusive 

7o Ln is an array of coefficients of the Lagrange polynomial 
Ldenom=l ; 
Lnom=l ; 
for i=l:N 
if i~=n 

Ldenom=Ldenom* (a(n) -a(i) ) ; 

Lnom=conv(Lnom, [1 ,-a(i)] ) ; 

end 

end 

Ln=Lnom/Ldenom ; 



function pl=polypower(p,n) ; 

7op: an array of coefficients of the Lagrange polynomial 
7. n: integer power 

7. pi: array of coefficients of polynomial p raised to the nth power 

pi=p; 

if n>l 

for i=l:n-l 

pl=conv(pl ,p) ; 

end 

end 



In order to verify the accuracy of the above implementation, we present the 
following arithmetic examples of the application of the above source code for 
the calculation of the semisimple part Sa, taken from [5]. 
" 3 4 3 " 
LetA= 2 7 4 
-4 8 3 

11 with multiplicity 2 and 1 respectively. The matrices Aq, Ai ac- 



It is easy to verify that A has 2 unique eigenvalues: 



ao = l,ai 
cording to Eq. 



are calculated as Aq 



1 

-0.2 



and Ai 



1 

0.2 



obtains the following values according, to Eq. 
The semisimple part is calculated according to Eq. 




1 



whereas the matrix c, 
1 11 
0.2 -2.2 
" 1 5.6 2.8 
Sa^ 8.6 3.8 
4.8 3.4 

calculating the Jordan decomposition of Sa which is a diagonal 3x3 matrix), 
as well as that {A — Sa) is nilpotent, since {A — Sa)"^ is the 3x3 zero matrix. 
The result obtained for Sa by the proposed method is identical to the result 



It is easy to demonstrate that Sa is semisimple (e.g. by 
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produced in [S]. 

The second arithmetic example is also taken from [5,. Let 



A 






-1 


























-1 



































1 


-1 




















0.25 





-1 


1 
































2 


-1 




















0.25 





-1 


2 
































3 


-1 




















0.25 





-1 


3 
































4 






















0.25 





-1 
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It is easy to show that the roots of the characteristic polynomial of A are the 
following: oq = 5,ai = 4, 02 = —1,13 = 3 (multiplicity equal to 2),a4 = 2 
(multiplicity equal to 2), 05 = 1 (multiplicity equal to 2),a6 — 0. The matrices 



Ao 



A3 



are calculated as: 




A, 



1 



A,; i = 0, . . . , 6 according to Eq. (|2.2 
1 1 . _ r 1 
-4.9 1 J' ^ ~ [ -2.5667 1 J ' [-1.667 

1 J ' ^ [ 1.1667 lj'^^[ 2.5667 1 
The matrix Cij obtains the following values according to Eq. 
-10 12 3 4 5 

-4.9 1.1667 -3.5 -10.2667 -24.5 
part Sa is calculated according to Eq. (|2.6|) : 



1 

4.9 



Ae 

The semisimple 






-1 


























-1 



































1 


-1 




















0.25 





-1 


1 




















-0.0156 


0.0156 








2 


-1 














-0.0156 


0.0156 


0.25 





-1 


2 














0.0039 


-0.0026 


-0.0156 


0.0156 








3 


-1 








0.0052 


-0.0039 


-0.0156 


0.0156 


0.25 





-1 


3 








-0.0004 


0.0002 


0.0039 


-0.0026 


0.0156 


0.0156 








4 


-1 


-0.0007 


0.0004 


0.0052 


0.0039 


0.0156 


0.0156 


0.25 





-1 


4 



Again this result is identical with Example 3 from [5] 
that the nilpotency holds, since {A — SaY — 0. 



It is also easy to verify 



5 Further generalizations of Theorem 3.1 and 
derived identities 

The main result of this section is the generalization of Theorem 13.11 which is 
expressed in the following Theorem. 
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Theorem 5.1. Let g G Ik[a;] be separable of degree n. Let H S (lk[a;])^"""" sue 
that det(n) ^ 0. Then for any polynomials fo, fi ,. . . , fm-i G there exist 
unique polynomial r £ k.[x] of degree less than mn, such that 



mod g, 





r 




fo 




r' 




fl 


n 




= E 






j.{m-l) 




fm — 1 



where E := gcd(g, det(n)). 
Proof. We start from 



nn = nn = det(n)/ 



(5.1) 



where 11 is the adjugate of 11 and /„i is the m x m unit matrix. 

Moreover, since E — gcd(g, det(n)) one has that there exist H,G <E k[x] 
such that 

Hdet(Il) +Gg ^ E. (5.2) 
Combining (|5.ip with (j5.2p we get, 



Hnn ={E- gG)I„ 



(5.3) 



Now, according to Theorem 13.11 there exist unique r 6 k[x] of degree less than 
mn such that 



r 




/o 


r' 




fl 




= HIi 




^(m-l) 




fm—l 



mod g. 



(5.4) 



Multiplying (|5.4I) from the left by 11 and setting (|5.3p in the resulting identity 
we get the conclusion. □ 



Remark 5.2. Choosing 11 = /„j m Theorem \5 . l\ we get Theorem \3.1\ Therefore 
Theorem \3.1\ can be regarded as a generalization of Theorem 15. il 



Now we will apply Theorem 15.11 to produce some formulas for polynomials 
involving derivatives. 

Corrolary 5.3. Let g, go, gi, ■ . ■ ,3m-i G k[x] be polynomials. Assume that g 
is separable and that 

{9i,9) = 1: < i < TO - 1. 

Then for any fo, /i, . . . , /m-i £ k[x] there exists unique r G k.[x] of degree less 
than mn, such that 



{i) - 



fi mod g, < i < m ^ \. 
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Proof. Applying Theorem 15 . 1 1 by 



n 



90 
91 







gr, 



and afterwards using that E = gcd(g, det(n) — Jlfc^o 5«) = 1 we directly get 
the conclusion. □ 

Corrolary 5.4. Let g, go,gi, . . . ,gm-i G ^[x] be as in Corollary \ 5. 31 Then for 
any /o, /i, . . . , /m_i S k[x] there exists unique r of degree less than mn such 
that 

(gjr)'^*' = fi mod g, < i < m - I. 



Proof Let 11 e (k[a;])'"^™ be the matrix defined by 



n 



50 







(5.5) 



/m-l\ (m-l) /m-l\ (m~2) 

I j.ym-1 I 1 ji/m-1 • • • \m-l)9m-l 

From the assumptions 

{9i,9) = 1: 0<j<m-l 

we have 

(det(n),.g)-l. (5.6) 

Applying Theorem 15.11 for 11, as given in (|5.5p . using (|5.6I) and the Leibnitz's 
rule, we get the conclusion. □ 
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